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It is shown that excitation spectra of Generalized Gibbs Ensembles (GGE) of one-dimensional 
integrable models with isotopic symmetry contain universal features insensitive to details of the 
distribution. Namely, the low energy limit of the subsystem of isotopic (for instance, spin) excitations 
is described by the effective action of a ferromagnet at thermodynamic equilibrium with a single 
temperature and with the stiffness determined by the initial conditions. The condition of universality 
is that the entropy per excited particle is small. 

I. GENERAL 

Dynamics of closed interacting many-body systems has lately became a subject of intense interest. Systems can be 
prepared in out-of-equilibrium state by any non-adiabatic change of the Hamiltonian (for instance, by quench) or by 
pumping. In such experiments many-body system is prepared in a state which is not an eigenstate of the Hamiltonian. 
It is then allowed to evolve coherently with many eigenstates evolving with different energies. In this paper we will 
focus exclusively on the fate of integrable systems. 

After some evolution time, which for integrable systems we expect to be of the order of the energy of the lowest 
energy excitation, the system can be considered as completely dephased. Then when calculating expectation value 
of any physical operator, it is reasonable to neglect off-diagonal terms between different eigenstates since such off- 
diagonal terms oscillate very rapidly in time and cancel each other. To be more rigorous one may consider averaging of 
expectation values over small time interval. Then it becomes possible to describe the system not as a pure state, but 
as a density matrix corresponding to some Generalized Gibbs Ensemble (GGE). Physical properties of the resulting 
GGE are determined by the integrals of motion which expectation values in the initial state fully characterize it 
[H , 01 1 Hi; 01 ; Hi- Conditions for validity of GGE for general integrable systems have been recently discussed in @. 

A priory there is no reason to expect any universality of resulting GGE. However, in the next sections we will demon- 
strate that GGEs appearing in quench experiments with one-dimensional integrable models with isotopic symmetry, 
have a universal structure that are insensitive to details of the distribution. Moreover, the low energy limit of the 
subsystem of isotopic (for instance, spin) excitations is corresponds to a ferromagnet at thermodynamic equilibrium 
with a single temperature. 

The defining feature of intergrable systems is existence of infinite set of mutually commuting operators Ij (integrals 
of motion). Alongside with the momentum operator Iq = P and the Hamiltonian I\ — H these conserving quantities 
can be extracted by expanding the logarithm of transfer matrix in spectral parameter 6: 

lnT(9)=P + i6H+(iO) 2 I 2 /2 + ... (1) 

Generalized Gibbs Ensemble (GGE) described above is given by: 

P = Z- 1 exp(-^/? J / J ), (2) 

j 

where Ij are integrals of motion and j3j are corresponding Lagrange multipliers which value is determined by the initial 
values of Ij . It maximizes entropy while taking into account the constraints imposed by the conserved quantities [l| , 
7]. This hypothesis has been proven in some particular cases 

In this paper we consider several typical examples of GGE in relativistic massive field theories possessing some 
isotopic symmetry. To work with such models is easier because their integrals of motion are known explicitly being 
fixed by the Lorentz symmetry and requirements of locality. However, as questions may emerge from the point of 
view of practicality of such models, we will first address prospects of their experimental realization. 



II. MODELS 



In this Section we give several examples of integrable models which find applications in cold atom physics. Some 
of them appear in the context of one-dimensional Bose gas. The situation relevant to our discussion is an experiment 
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where a one-dimensional trap containing condensate is abruptly split to become a double- well potential as, for instance, 
was done in Q . We will be interested in the situation when the split is not complete such that there is still a substantial 
inter-well tunneling. In the low energy limit interacting ID Bose condensate is uniquely characterized by its local 
phase field <p(t, x) which dynamics is described by the Tomonaga-Luttinger liquid model (the model of non-interacting 
bosonic field). When the trap is split there are two condensates and an additional term appears in the Hamiltonian 
which describes their coupling. As has been demonstrated (see e.g. [HI), the Hamiltonian describing a system of 
two parallel traps consists of two independent parts. One describes the symmetric part of the condensate phase field 
<j) + = (4>i + <j)2)/y/2 and the other one describes the asymmetric part <f>_ = (</>i — (j> 2 )/y/2: 



H = H 



L 
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d.x 
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K-\d x 
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K(d x <j)-) 2 +K- 1 {d x 6-f - Acos(V27rc/>_) 



(3) 
(4) 



where v is the phase velocity, K is the Luttinger parameter normalized in such a way that in Tonks-Girardeau limit 
K = 1, and parameter A is proportional to the inter-trap tunneling. Fields 4> a ,6a obey the standard commutation 
relations: 



[4> a {x), 9 b (y)] = S ab Q H (x - y), 
where Qh (x) is the Heaviside function. The low energy limit of the original bosonic fields is given by 



1>i 



(5) 



(0) 



where dots stand for operators with higher dimensions. 

Models (I3l4p are integrable; the first one describes a free field, the second one is the famous sine-Gordon model. 
For K > 1/4 the cosine term in ((4]) is relevant and scales to strong coupling. As a result the sine-Gordon model has 
massive spectrum of the relativistic form: 



E n (p) = y/(vp)*+M*. 



(7) 



Its excitations include solitons and antisolitons corresponding to kinks interpolating between neighboring minima of 
the cosine potential and (for K > 1/2) their bound states called breathers. Their masses in term of the soliton mass 
Mo are 



M n = 2M Q sin 



2{AK - 1) 



l,2,..[4if - 1] 



(8) 



Sine-Gordon model (j4} with an additional chemical potential field also describes the superfiuid-insulator transition (llj. 
This fact is important for practical applications since in cold atom physics there are well developed techniques for 
observation of such transition and of its associated features [l2l|. 

Another relevant model of ultracold atom physics is the model of SU(N) fermions with a point- like attractive 
interaction. The corresponding Hamiltonian is 



(9) 



where a,b = 1,2, ...N. Experimentally one dimensional systems of SU(2) fermions have been realized by Moritz et 
al. [13 ■ Fermions with SU(N) symmetry have been recently realized in optical lattices in [l4|, IHj]. Earlier theoretical 
work on SU(N) fermions focused on their equilibrium properties [l6l - fl9| and their applications to quantum information 
processing|20|. 

In the limit of weak coupling g << (fi/m) 1 / 2 model ([§]) becomes equivalent to another famous integrable model, 
the the so-called Chiral Gross-Neveu one. In this limit one can linearize the spectrum of fermions near the Fermi 
points and replace the operators: 



tP(x) = e- lkFX R(x) 



,lk F x 



L(x), 



(10) 



where slow fields R,L contains Fourier harmonics with momenta much smaller than the Fermi momentum hp. Sub- 
stituting pU|) into ([§]) one obtains 



H = J dx(^- ivR+d x R a + ivL+d x L a - gR+L a L^R^ , 



(11) 
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where v — kp/m. 

Both models © and its relativistic limit (JTTJ) are integrable and together with the sine-Gordon model (UJ) are 
among the best studied models of that kind. At g > model ([TT|) is asymptotically free with the interaction scaling to 
strong coupling in the infrared limit. The spectrum of (jlip is split into two independent parts, one of which remains 
gapless and will not be discussed. The spin sector having SU(N) symmetry have massive excitations with spectrum 
(0 consisting of the fundamental particle with mass Mq ~ ^(g/v) 1 ^ exp(—2irv/Ng) and its bound states 

M n =M «, n = l )2 ,..JV, (12) 
sin(7r/7V) 

These excitations all carry isotopic indices of the SU(N) group and transform according to its fundamental represen- 
tations described by vertical column of Yang tableau of length n. 

Model is particularly interesting for us since it has a non-Abelian continuous symmetry which plays a principal 
role in the subsequent discussion. 

III. TBA EQUATIONS FOR GGE 

Models (|4llip were among the first field theories solved by the Bethe Ansatz (BA). To simplify our consideration 
we discuss in detail only the simplest case of SU(2) invariant Gross-Neveu (GN) model N = 2 and also touch on the 
sine- Gordon model (|3]). 

A. General facts about Bethe ansatz 

Relativistic (Lorentz invariant) models are particularly convenient for our discussion since they admit a simple 
classification of integrals of motion. This classification becomes particularly transparent in the so-called rapidity 
representation when energy and momentum of particle of mass M are parametrized as 

E = McoshO, P = Msinh9, E 2 - P 2 = M 2 , (13) 

(we set v = 1 for simplicity). Then Lorentz transformation (boost) becomes just a shift of rapidities of all particles: 
6i — > 6i + a. Consequently, all integrals of motion can be classified according to their Lorentz spin: 

£7±P = m5> ± * lH\±ljl\=M^ Ml+m , J = 1,2,... (14) 
i=i i=i 

Being a Lorentz invariant object the A-body scattering matrix depends only on difference of rapidities of individual 
particles. For integrable models such S-matrix can be written as a product of two-body scattering matrices. Until 
the system is somehow restricted (for instance, placed in a box), particles rapidities in (|14l) are arbitrary. However, 
as soon as motion of the particles is restricted, this changes. 

The Bethe ansatz (BA) equations which determine eigenvalues of all integrals of motion for a model in a box of 
length L with periodic boundary conditions can be derived from the straightforward solution [2l[ or taking the SU(N)- 
invariant solution of the Yang-Baxter equations for the two-particle scattering matrix and applying the methods of 
factorized scattering [22j . The condition for the periodicity of the wave function of relativistic interacting particles of 
mass M is 

M 

exp(iMZsinh0 i )£ = - 0j)£, (15) 

where £ is a vector depending of spin indices of the particles. The meaning of this equation is straightforward: the z-th 
particle going around the system scatters on all others (it does it one-by-one which is the condition of integrability) 
acquiring a phase factor given by the product of all ^-matrices on the right hand side of (|15[) . This phase factor is 
compensated by exp(ipiL) in the left hand side of this equation. 

For models with internal (isotopic) symmetry S is a tensor and diagonalization of (|15|) requires some effort. The 
result of this diagonalization is the so-called nested Bethe ansatz. For the GN and the sine-Gordon models the result 
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N 



M 



exp(iAfisinh^) = ]JSo(9 t - Oj) JJ e x (0< - A Q ) 



M 



l[e 1 (X a -e i ) = Y[e 2 (X a - A h ) 



(16) 
(17) 



»=i 



6=1 



where for the GN model 



and for the sine-Gordon it is 



e n (x) 



x — imr/2 
x + in7r / 2 



e n (x) 



sinh 




- irar/2) 


sinh 




+ iwr/2) 



7= (4if- 1)" 



where Sq(9) is some known function which exact form is not important for the present discussion. Number J\f stands 
for the number of particles and M. is related to the spin projection: S z = Af/2 — A4. The qualitative difference 
between the GN and the sine-Gordon model is that in the latter case for K > 1/2 function So (9) have poles on 
the physical strip and there are bound states (breathers). These breathers carry no spin. We have to add that the 
integrals (|14|) with integer I are local (they have integer Lorenz spin /); in principle, one can imagine that quench 
generates non-local integrals with non-integer I. However, the further discussion does not depend on whether it is 
true or not. 

Generalization of the nested BA equations (116I17[) for models with other simple Lie group symmetry follows the 
standard scheme described, for instance, in [23|]. For a given simple Lie group one has to modify Sq(9) and replace 
(jTTJ) with a hierarchy of coupled algebraic equations for rapidities A^ (j =1,...N-1) where ^ — 1 is the dimension of 
the corresponding Kartan subalgebra. The structure of the hierarchy reflects the structure of the Dynkin diagram for 
the given group. 



B. Derivation of the universal dynamics 

The emergence of the universal spin dynamics can be ascertained already from (|16I17I) . Since the particles are 
massive, their number in the ground state is zero: Af = M. = 0. However, a nontrivial GGE emerges after a work had 
been performed on the system resulting in a finite particle density. For the infinite system the values of all integrals 
of motion are determined by the rapidities 9i\ the auxiliary variables A appear only when the system is put in a box. 
Therefore, at least in the limit of small particle density, one can neglect a feedback on 0's from A's and consider the 
distribution of (9's as an independent function. 

When the ratios Ij / L are finite in thermodynamic limit the distribution function of rapidities 9 must decay suffi- 
ciently fast at infinity. Then in the limit of large A's one can replace in Eq. (fl7|) 

M 

nei(A„-0i)« [ei(X a )f (18) 

i=l 

which indicates that the spin sector (described by A's) decouples from 9's. In fact Eq. flTTj) with substitution (fT8)l 
resembles the BA equation for a spin S=l/2 Heisenberg magnet. In order to figure out whether this is ferro- or 
antiferromagnet more detailed analysis is needed (see below). As we will show, it is ferromagnet. The difference 
between the SU(2) GN model and the sine-Gordon model is evident already at this stage: in the former case we have 
an isotropic magnet and in the latter case it is anisotropic (U(l) or easy plane magnet). 

In order to see the emergence of universal spin dynamics and establish its conditions we have to derive Ther- 
modynamic Bethe Ansatz (TBA) equations. This derivation follows the standard scheme. First, we establish that 
generically complex solutions of Eqs. pT|) in the thermodynamic limit (L —> oo, M/L, M./L = finite) have only 
fixed imaginary parts. More specifically, these solutions group into clusters with a common real part (the so-called 
"strings"): 



X n ,j:a = f£ n) +m(n + 1 - 2j)/2 + 0(exp(-constL)), n = 1,2, j = 1,2, ...n. 



(19) 
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Then we introduce distribution functions of rapidities of string centers p n {9) and the distribution function of particle 
rapidities po(0). Functions p n (9), po(9) describe distribution of unoccupied spaces. Their ratios are parametrized by 
excitation energy functions e„: 

p n {9)/p n {9) = exp[-e n (0)], p {9)/ Pa {9) = exp[-e o (0)], . (20) 

/CO p 
d6 Po (6), M/L = Y,n d6 Pn (6). (21) 
n=l •* 

The entropy of the state is given by the expression 

oo . 

S = L / d # [(fin + Pn) ln(p„ +Pn) — Pn I» Pn ~ Pn In Pn ] ■ (22) 

71=0 ^ 

The TBA equations are a result of minimization of the generalized free energy 

fl = Y^Pjlj -S = J d9K(9) Po (9) - S (23) 

3 

where 

K{0) = Y,\Pke M +P k e- ke ], (24) 
fe 

subject to constraints imposed by the equations for the distribution functions: 

M 

Pn + Pn = S * {pn-l + Pn+l) + ^.OT" COSh6», (25) 

Z7T 



The result is 



c 



n {9) = s * ln[l + e^ e ~>][l + e £ »+'W] - S n>0 K(9) (26) 
d9_ 
2^ J 



n/L = - I —K{9) In [l + e £oW l (27) 



s * f(9) 



A9'f{9') 



, 2?rcosh(0 - 9') 
From these equations one can restore the integrals of motion: 

A peculiar property of TBA equations (|25I26[) is their quasi- locality: given e n ,p n are related only to their nearest 
neighbors. Therefore, if €q,Po are fixed the TBA for e n ,n = 1,2, ... will have the same form as (|2"6")l . but with K 
replaced by G{9) = s * ln(l + e eo ^). 

It is natural to assume that all integrals of motion are extensive quantities ~ L and hence their densities are finite. 
It follows then from (|27|) that eo(0) — > — oo at infinity (changing sign more than ones at finite 9). Then the function 
G{9) is integrable which will be important for what follows. On Fig. 1 we give an example of such function. We 
have chosen eo in such a way, that G{9) has some nontrivial structure at small 9. This structure is determined by the 
initial conditions of the quench. All these details, however, does not make a difference in in the asymptotic region of 
large 9 (see Fig. 2) where the bahavior of e„ is determined by a single integral characteristic of the ^-distribution (see 
the derivation below). 

Inverting the kernels in (I25I26[) with n > 1 we obtain TBA equations for GGE in the following form: 

ln[l + e 6 »W]- A„ m *ln[l+c- e '"W] =a n *G{9), n, to =1,2,... (29) 

Pn + Aim * Prn = »n * S * p (30) 

where the Fourier images of the kernels are 

a n (u>) = cxp[— 7rn|cj|/2], A nm (uj) — coth(7ra;/2) exp(— \n — m|7r|a;|/2) — exp(— (n + m)Tr\u>\/2) . 
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FIG. 1: An example of integrable function G(9; {/?}). Function eo changes sign six times going to — oo at both infinities. 




FIG. 2: e„(9) functions for n — 1, 2, 3 corresponding to G(8) shown at Fig.l. The asymptotics does not depend on fine details 
of behavior G(9) at finite 9. 



Eqs. (|29l30[) with G and s * po replaced by the delta functions are precisely TBA's for a spin S=l/2 Heisenberg 
ferromagnet in thermodynamic equilibrium (!) (24|. When these functions are not delta functions, but just sharp 
peaks, the analogy with the ferromagnet remains valid just asymptotically. For this analogu to hold it is sufficient 
that the integral J G(y)dy converges, but is >> 1. This would correspond to low effective temperature limit of 
the fcrromangnet when the free energy is determined by large rapidities (small momenta) such that details of the 
dispersion at large momenta are not important. Indeed we have 

J/T=±JdyG(y). (31) 

A similar condition must be satisfied by the function po(0): 

J dyPo(y) = finite. (32) 

The condition J G(y)dy » 1 means that there is little feedback from e„, n > 1 to eo & n d therefore one may use 
function G(9) to characterize the state. This is more convenient than to use the intergrals of motion which themselves 
can be restored from G via |28|) . From (f27|) and (j3"T|) it follows that the effective exchange integral of our ferromagnet 
is of the order of the energy density per particle. Since the entropy per particle in the ferromagnet (T/J) 1 / 2 , the 
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requirement T/J « 1 is equivalent to requiring the entropy per particle to be small. Under that condition the 
universal spin dynamics emerges in GGE which is, perhaps, the most striking result of our derivation. It is also clear 
that the spin subsystem of GGE is at thermal equilibrium and is described by the ordinary Gibbs ensemble with a 
single temperature. 



IV. CORRELATION FUNCTIONS 



To understand how the universal dynamics described in the previous subsection is related to observable quantities 
one has to consider correlation functions. This is a difficult problem and though we are not in a position to offer a 
detailed solution, we feel obliged to make some remarks. 

To be closer to real experimental systems we consider correlators of the bosonic creation and annihilation operators 
in the model of two coupled Bose condensates (I3l4j) . In the low energy limit these operators are expressed in terms 
of the phase fields <j>± ([6]). Since the two sectors of the model are decoupled, the correlation functions factorize. For 
instance, for the 2-point one we have: 

(^ a (M)^ + (0,0)» ~ po (( c iV^^(*^)e- i V^* + (o.o)))(( c i ^^(*^) e ( 1 - 25 -) i ^'A-(o^) )) (33) 

The correlation functions of bosonic exponents of gapless fields in GGE have been calculated [25| , but correlators of 
the sine-Gordon fields have never been analyzed in this context. Below we will discuss some general features of 
these correlators. 

We will discuss the limit of zero effective temperature T e ff = when all e n — > oo in (|29p . but their ratios remain 
constant. This is the limit when GGE distribution is reduced to a single "vacuum" state. This state is characterized 
by some distribution of the particle rapidities and is ferromagnetic, that is has maximal possible spin. Using the 
Lehmann expansion where Green's functions are expanded in matrix elements over excited states \n): 

(0(t,x)0 + (0,0)) =J2\ <n|O(0 J Q)|«oc> [V"^-®-" )-^- p -" ) (34) 

n 

In the sine-Gordon model excitations are classified as solitons, anti-solitons and their bound states (see Section II). 
(Anti)solitons can be described as particles with (negative) positive spin projection (s — ±1/2). The ferromagnetic 
state in this context corresponds to the state with only one type of particles (for example, solitons). We denote such 
state as 

l(*n,.*) s ;(o) AS > 

According to '2(1 . the following matrix elements do not vanish: 

< (0)^|e^- ] .-^ +JV ) s ; (^x, -^)^ >= 

<vac\e 1 ^-\(9' 1 ,...e' n+N ) ; (Qx-Jn^ + U ...m + 6 n ) > (35) 

V / $ \ / AS 

These matrix elements would correspond to excitation of N magnons and therefore the spectral function of bosons 
contains multi-magnon processes. Unfortunately, the form is these matrix elements is very cumbersome which makes 
further calculations difficult. At the moment the moment this is all we can say. 



V. OTHER MODELS 



The results obtained for the SU(2)-invariant model (jTTJ) can be easily generalized for any simple Lie group. Qual- 
itative differences appear only when the S-matrix for physical particles contains the so-called RSOS (Restricted 
Solid-on-Solid) component. Such models have an exotic degenerate ground state and excitations with non-Abelian 
statistics. In view of rarity of such problems we will confine ourselves to a brief discussion. 

A typical representative of this class of models is the SUfe(2) Wess-Zumino-Novikov-Witten (WZNW) model per- 
turbed by the current-current interaction term: 



H 



d.r 



2ttv 



j a r 



j a j a 



gJ a J° 



(36) 
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where the current operators satisfy SU(2) Kac-Moody algebra of level k: 

[J a (x), J b (y)} = ie abc J c (x)S(x - y) + ^5 ab 5'(x - y) (37) 

Z7T 

which, in fact, coincides with the commutation relations of the fermionic bilinears 

k k 

J a = ^ R j~a a aP R 3P> J" = X] L ta a aff L j0 
3=1 J = l 

Model (|3^|) describes the SU(2)-invariant sector of the SU(2)xSU(k) model of fermions [27j ■ The soliton excitations of 
this model carry zero modes of parafermions; this is the origin of the non-Abelian statistics (see, for instance [28], [29] 
for the discussion). This model also emerges as a continuum limit of the lattice model of spin S = k/2 integrable 
magnet with a small Ising-like anisotropy A[30j|: 

H = J2 Vk [i S t S I+i + h.c.),S?S] +1 ;A] (38) 

3 

where 



P fc (x,y;A)= X A nm (A) 



x n y m 



n+m<k 



is some known polynomial. The coupling constant in (|36p g ~ y/ A — 1. 

The TBA equations for this model differ from (|26I24I25[) only in one respect: the driving term is placed not in the 
first, but in the fc-th equation and the free energy is also determined by et.[27j. Therefore inverting the kernels in 
TBA yields not one, but two sets of independent equations: 

ln[l + e £ "+ fc(0) ]-A nm *ln[l + e- £m + fe(e) ] =a n *G(9), n, to =1,2,... (39) 
ln[l + e £ " (e) ] - A nm * ln[l + e~ em ^] = B n * G{8), n, m = 1, 2, ...k - 1 (40) 



where 



sinh[7r(fc — n)oj/2] sinh[7r(fc — max(n, to))oj/2] sinh[7rmin(n, m)w/2] 

o„ = . 7-. , — 7777 — , A nm = 2coth(7ro;/2)- 



sinh(7rA;a;/2) ' sinh(7rfca;/2) 
and the corresponding equations for the densities are 

Pn+k + A nm * p m +k = «n * S* p~ k (41) 
Pn + Anm * Pm+1 = B n * S * pk (42) 

The set of equations for n > k describes spin S=l/2 ferromagnet as before. As far as the equations for n < k are 
concerned, in the low energy limit they describe a conformal theory. From (|40[) and (|4"2"]) it follows that at large |0| 

e n (6)~e-^/ N , ~p n {e)~e-WN^ (43) 

and, since p(8) = 2ir f A6' p(6'), the spectrum is linear. The low energy limit of the second one corresponds to the 
conformal field theory of parafermions. 



VI. PHYSICAL CONSEQUENCIES 

We conclude this paper by reiterating its main result: provided the initial state of our system containing a gas of 
excited particles is a low entropy state with (small entropy per particle) , its spin dynamics is universal. This condition 
is equivalent to the condition T << J (the effective exchange integral ([3T10 . Namely, in the low energy limit GGE 
spin subsystem looks like a ferromagnet at thermodynamic equilibrium. GGE's of models with non-Abelian statistics 
also contain a critical sector described by some conformal field theory (for the example given in the main text it was 
the theory of Zj, parafermions). A qualitative explanation for this result is simple: quench creates particles not fixing 
the total magnetic moment and the system chooses a state with maximal entropy which is a state with maximal total 
moment. 
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There two natural questions to ask. The first is how to prepare such a state, the second is how this universal 
dynamics reveals itself in observable quantities. One possible answer to the first question is to use pumping with an 
appropriate field generating pairs of excitations in the spin sector. In the case of ultracold atoms in optical lattices 
such pumping can be performed with optical lattice modulation. A pulse with duration to and frequency uj may excite 
pairs of particles of type n momenta provided that ui > 2M n . Assuming that pumping does not give momentum to 
the system, excitations will be generated as pairs of quasiparticles with opposite momenta. Energy conservation then 
gives that such quasiparticles will be excited with rapidities centered around M n cosh 9 ss w and the width of energy 
distribution ijj" . 

Now let us discuss the second question. What are the physical consequencies of the effectively "thermal" character 
of the spin subsystem? Obviously, one cannot argue that the form of generic isospin correlation functions will be 
that of a corresponding ferromagnet. This is clear from the analysis of Section IV. Therefore even though we find 
the same distribution of eigenstates as in the appropriate ferromagnet, matrix elements of the operators may be 
very different. Is it possible then to perform any physical experiments, which would demonstrate the corresponding 
" thermal" character of the isospin sector? 

The fist option is to measure fluctuations of the magnetization. For the Chiral Gross-Neveu model the former 
corresponds to selecting any of the SU(N) operators and measuring the value of this operator in a finite segment of 
the svstem pjll 132], Unlike more generic operators, smooth components of the spin density operators have the same 
matrix elements for the Chiral Gross-Neveu and ferromagnetic models (analogously A_ ~ d x <j>- operator for the 
sine-Gordon has the same matrix elements as the S z operator for the easy-plane ferromagnet). We emphasize that 
measurements of the spin operator should be repeated many times so that one could extract not just the average value 
but all fluctuations of the operator. Putting it differently one can say that results of individual measurements should 
be combined into distribution function. We predict that the distribution function of magnetization fluctuations of the 
Chiral Gross-Neveu model after a quench should be the same as in a ferromagnet at finite temperature. 

The second way of observing thermal character of the isospin sector is to measure its fluctuations of energy. In a 
thermal ensemble energy fluctuations are given by the specific heat (AE 2 ) = cyT[33|. In equilibrium specific heat 
is itself a function of temperature. Hence by measuring the average energy and its fluctuations one can effectively 
measure the equation of state. Our analysis suggests that the average energy and fluctuations of the energy in the 
Chiral Gross-Neveu model following the quench should be given by the equilibrium equation of state of an SU(N) 
ferromagnet. Measurements of energy fluctuations in the system (or in a fragment of the full system) can be done 
experimentally (see e.g. experiments [34| . In these experiments the average energy of an interactin g ID Bose gas 
was measured in an array of tubes. Recent experiments allow local resolution of individual Id systems [35H37j ] . which 
should make it possible to measure not only the average energy but also energy fluctuations. Assuming separation 
of the isospin and density sectors, it should also be possible so separate the isospin part of the energy from the total 
energy. 
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